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A generalized Green-Kubo formula is derived for a quantum dissipative system of driven 
Brownian particle, in which the coupling between the system and the environment is linear. 
The structure is essentially the same as that for the generalized Green-Kubo formula for 
driven granular particles. It is demonstrated that the correction to the conventional Green- 
Kubo formula is zero for a free Brownian particle. 

§1. Introduction 

Green-Kubo formuldD' is one of the most fundamental relations in nonequilib- 
rium statistical physics. The original derivation is restricted to the linear nonequilib- 
rium case, but number of generalizations are proposed by many researchers P''^'®''^''^ 
though the relationship among their formulations is not still well understood. 

We believe that roles of dissipation in Green-Kubo formula should be clarified, 
though it is unclear in the original derivation. Indeed, the Green-Kubo formula 
defines the transport coefficient which represents the dissipation. Moreover, if there 
is no dissipation in a system, the time integral of current correlation function in 
the formula should diverge. Therefore, purely mechanical derivation of Green-Kubo 
formula might be misleading, but correct derivation should include the dissipation 
explicitly. 

Recently, Chong et al.*^ have derived a new generalized Green-Kubo formula 
for driven dissipative and classical particles as a natural extension of that by Evans 
and Morriss.'^ Their derivation has several remarkable points; (i) the formulation can 
include the integral fluctuation theorem without microscopic time-reversal symmetry, 
(ii) the role of dissipation is clear in their derivation, and (iii) one can develop the 
nonequilibrium mode-coupling theory for sheared granular liquids or sheared glassy 
systems. Since their derivation is so general that one can expect that their method 
can be used for quantum cases. In this paper we demonstrate how to apply such a 
formula to quantum systems. 

In this paper, we focus on a nonequilibrium steady state of a quantum Brow- 
nian partilce. This is because (i) the simplest system among open quantum sys- 
iem.s^^^^^^^^ (ii) we know the origin of dissipation of quatum Brownian system 
as the energy flux between the system and the environment, and (iii) there is the 
quantum version of the violation of fluctuation-dissipation relation.^I^ 

The organization of this paper is as follows. In the next section, we specify the 
basic equations to be analyzed in this paper. In section 3, we will obtain an exact 
solution of our model. In section 4, we will present the formal form of generalized 
Green-Kubo formula whose bilinear form is reduced to the conventional Green-Kubo 
formula. In section 5, we will focus on the case that a Brownian particle in a harmonic 
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potential, in which the integration involved in the generalized Green-Kubo formula 
can be carried out. In section 6, we will discuss and conclude our results. We also 
invole two Appendices, where Appendix A is devoted to the derivation of generalized 
Kubo's identity, and Appendix B gives an explicit calculation of the time evolution 
of the momentum of the Brownian particle. 



§2. Model 



Let us begin with a quantum master equation for a Brownian particle^' under 
an external field F^x, which is the essentially same as Caldeira-Leggett model.'^ 
Note that the generalization to N Brownian particles is straightforward from that 
presented here. For simplicity, we restrict our interest to a particle coupled with the 
heat bath. 

Let us consider a Brownian particle of mass m with its coordinate x and mo- 
mentum p in a potential V{x) under a steady external field -FexEE The particle is 
assumed to be coupled with a bath consisting of a large number of harmonic oscil- 
lators with masses m„ and frequencies tUn- Thus, the total Hamiltonian might be 
written as 

H = Hs + He + Hb + Hi + Hex 
Here, the system Hamiltonian of the particle is given by 

2 

Hs = ^ + Vix), (2-2) 

where we do not specify the form of the potential in the formulation in the main 
part of our paper. The bath Hamiltonian is represented by 

HB = ^fuJn Uibn + =Y1 (2^^" ^'r^nUjlxlj . (2-3) 

Here, bn and bn denote the Bosonic creation and the annihilation operators of the 
bath, respectively, while x„ and pn are the corresponding coordinate and the mo- 
mentum. Similarly, the interaction Hamiltonian is given by 

Hj = —X KnXn = —xB, (2-4) 
n 

where Kn is the coupling constant, the bath operator is 



B = y~^ HnXn = Hn\l 7r~ (^n + (2-5) 



*-* Of course, the external field can be involved in the potential, but we separate the contribution 
of the external field from the stationary potential V{x). 
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We introduce the external Hamiltonian coupled with the external force -Fex as 

Fex = -xFe^. (2-6) 

Note that //gx can be absorbed in Hs but we separate its contribution to clarify the 
response to the external field. We also introduce the counter-term Hamiltonian: 

This counter-term Hamiltonian can be absorbed in the potential term V{x) as 
Ve{i{x) = V{x)+Hc- However, if we regard the unperturbed Hamiltonian as Hs+Hb 
without couping between the system and the bath at t = 0, the effect of He appears 
in later expressions. We also note that He must be treated as a term second order 
in the coupling, while Hj is of the first order. 

Caldeira and Leggetl!!^ were interested in the low frequency bahavior of parti- 
cles. Then they adopted the simple assumption Hujq <C Mm{hil, InkT} under the 
condition V{x) ~ (l/2)ma;QX^ + O(x^) with a high frequency cutoff i7, the Boltz- 
mann constant k and the temperature T: 

xs{-t) = e-'^sr/h^^iHsr/h ^ ^ _ hHs,x]T. (2-8) 

h 



Here, however, we do not have to use this quasi-classical expression (l2-8p for our 
argument. Here, we have introduced the commutation relation [A, B] = AB — BA. 

The starting equation is the Born-Markov approximation for the reduced density 
matrix of the Brownian particle, which obeys^' 



i iF 1 f°° 

;Ps{t) = --[Hs+Hc, ps{t)]+^[x, ps{t)]-- dTtTB[Hi,[Hi{-T), ps{t)0pB]]. 



d_ 

dt' 

(2-9) 

Hereafter, we adopt the interaction picture with the respect to the unperturbed 
Hamiltonian Hq = Hs + Hb ■ We shall assume that the initial condition satisfies 



PS{0) = Peq® PB, (2-10) 

where 

exp[-PHs] exp[-(3HB] , . 

with fi = 1/kT. Note that the assumption on ps{0) might be removable. Indeed, the 
steady distribution ps{t — > oo) of Caldeira-Leggett model with Eq. (j2-8p is relaxed 
to ps{0)'^'^ if the particle is trapped in a potential. 

For the discussion of quantum Brownian motion we introduce the spectral func- 
tions 

D{T)^i{[B,BB{-T)])B, D,{t)^{{B,Bb{-t)})b, (2-12) 
where ^^(t) = e'^Bt/h^^-iHBt/h {A, B} = AB + BA. Note that Hs and Hb 
are decoupled with each other in the unperturbed state, where {■■■)b represents 
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the average in terms of the density matrix ps- We also note that D{t) and -Di(t) 
are respectively referred to the dissipation and the noise kernel. Making use of the 
spectral density 

J{co)=y2,r^S{u;-u;n), (2-13) 



we can write the explicit representations for the correlation functions 

/>oo 

D{T) = 2h dio J (uj) sin LOT, (2-14) 
Jo 

Di{T) = 2hj^ dujJ{uj)coth(J^^ cos ujT. (2-15) 
After straightforward calculation, eq. (|2-9p can be rewritten a^ 

^^psit) = -l[Hs + Hc,ps{t)] + '-^[x,ps{t)] 

+^ /, dr{iD{T)[x,{xs{-T),ps{t)}]-Di{T)[x,[xs{-T),psm)i2-lQ) 

which is the basic equation of this paper. 

The properties of the second line of (I2-16P strongly depend on the behaviour of 
the dissipation and the noise which are determined by J{uo). We adopt a continous 
distribution of the bath modes and replace the spectral density by a smooth function 
of u; for the explicit calculation. 

For the calculation in section 5, we should specify the form of J {to). Here we 
adopt the Ohmic spectral density with Lorentz-Drude cutoff function, J{uj) is given 
by 

J ^ = '-u 2-17 

where 7 is a damping constant and 17 is a high frequency cutoff. Note that J{uj) 
satisfies Ohmic dispersion J{uj) — > 2m7u;/7r as w ^ 0. In this case, D{t) and Di{t) 
are respectively given by 

D(t) = 2m-fhn^e-^^ (2-18) 

and 

D,{r) = Am^kTQ^ J2 ^ J (2-19) 

n=~oo " 

for r > 0, where Un = 27mkT/h is known as the Matsubara frequency.'^ 
§3. The solution of Liouville equation 
Let us rewrite Eq. ()2-16p as 

f^psit) = -iC^psit), (3-1) 
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where the Liouville operator is given by 

i iF 1 1'°° 

iC^^-[Hs+Hc,]-^[x,]-^ dT{iD{T)[x,{xs{-T),}]-D,{T)[x,[xs[-rm. 

(3-2) 

Since this Liouville operator iC"^ is independent of time, we can use the identity 

^-^ch = 1 + / dse-"^''\-id^). (3-3) 

^0 

Substituting Eq. (|2-10p into Eq. (j3-3p we obtain 

ps{t) = peq + r dse-"^'\-iC^)p,^. (3-4) 

Here, —iC^peq consists of four terms: 

iF i i f°° 

- i/:Veq = -^[a::,/Oeq] - ■^^[x^,Peq] + ^ / dr D{t)[x , {x s{-t) , Pc^^]] 

1 

~2h? Jq dTDl{T)[x,[xs{-T),Peq]]. (3-5) 

The first term on the right hand side of eq. ()3-5p produces the conventional Green- 
Kubo formula. With the aid of Kubo's identity (|ATp . we can rewrite 

[X, peq] = Peq / dXc^""^ [Hs, xje'^^^ 

Jo 

ih ff^ 

= Poq / dXps{-ifiX), (3-6) 

Jo 

where we have used [Hs,x] = —ihp/m. Thus, the first term on the right hand side 
of (13-51) is reduced to 







^X,Peq\ = Peq / dXps{-ihX). (3-7) 

Similarly, the second term on the right hand side of (j3-5p which is the order of 
square of the coupling constant can be calculated. From ()A-ip we readily obtain 



[x\peq] = dX{xs{-ihX),ps{-ihX)}, (3-8) 

where we have used [Hs,x'^] = —ih{xp + px)/m. Thus, the second term on the right 
hand side of (|3-5p is reduced to 

- '-lC[x\p,^] = r dX{xs{-ihX),ps{-ihX)}. (3-9) 
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The contributions from the third term and the fourth term on the right hand side 
of (|3-5p are more comphcated. From ()A-ip . (\A-2\i and (|A-3P we obtain an identity 

[x,{xsi-T),Pcq}] = Pcci{[x,Xs{-t)] + [x, Xs{-ihP - t)]) 

-Peq r dXpsi-ihX){xs{-r) + xsi-ihp - t)). (3-10) 



m 

Similarly, thanks to ()A-ip we obtain 

[X, [Xs{-T),peq]] = Peq {[x,Xs{-ihP - t)] - [x,Xs{-t)]) 

r/3 



m 



ih^ I d\ps{-ihX){xs{-ih(3-T)-xs{-T)). (3-11) 







Thus, the third term and the fourth term on the right hand side of (j3-5p is reduced 
to 

^ Jo ~ 2^_/g dTD{T)[x,[xs{-T),Peq\] 

= r dTD+{T){[x,Xs{-T)]-- f dXps{-thX)Xs{-T)} 

- r dTD^{T){[x,xsi-f)] - - / dXps{-thX)xs{-f)}, (3-12) 
2ft^ Jo m Jq 

where f = r + ihp and 

D±{T) = Di{T)±iDiT). (3-13) 

Thus, we can write 

Ps{t)=Pe^+ dse-'^''[pec,0], (3-14) 
Jo 

where 

= ^ dXpsi-ihX) -— r dX{xs(-ihX),ps(-ihX)} 
m Jo m Jq 

/ dT{b-{T) / dXps{-ihX)xs{-T - ihf3) 







2hm 

D+{t) / dXps{-ihX)xs{-T)} 



+ 2^/ t^T(£'+(r)[x,a;5(-r)]-D_(r)[x,X5(-r-i;i/3)]). (3-15) 

Before closing this section, we should note an important property of which 
satisfies 

(0)eq = tr5{pcq0} = 0. (3-16) 



This relation is easily verified from Eq. (j3-5p with the invariant property of the trace 
under a cyclic permutation. 



Generalized Green-Kubo formula for a dissipative quantum system 



7 



§4. Generalized Green-Kubo formula 

Let us derive the generalized Green-Kubo formula. For simplicity, we discuss 
the average behavior of the momentum p 

{p)t = tvs{ps{t)p}. (4-1) 
Substituting ()4-ip into (j3-14p we obtain 

{p)t = trs{pcqP} + /* dstis{e-'^''[pcc,e]p}. (4-2) 
Jo 

With the help of the property 

tvs{e-'^''[pe^O]p} = tisiPc^OpHit)} (4-3) 

with pnit) = e'^^^p = g^i-^i/^p^-^^i/f^^ we can write the generalized Green-Kubo 
formula, where the Liouville operator iC satisfies Heisenberg's equation of motion 
^H{t) = iJ-^Anit) = jr^[H,AH{t)]. We note that the last equality for pnit) holds 
because of our special set-up, where we adopt the basic model under the Born- 
Markovian approximation with the linear couping between the system and the envi- 
ronment. 

The explicit form of iC is given by 

iC = ^[Hs + Hc,]-'^[x,] 

r {iD{T){xs{-T), [x, ]} + D,{r)[xs{-T), [x,]]) , (4-4) 



where we have used the invariant property of the trace under a cyclic permutation. 
We can itroduce 

A = iC^ - iC, (4-5) 

whose expectation value corresponds to the phase volume contraction in classical 
situations. The operator A is immediately obtained as 

I dTD^{T){[x,[xs{-T),]] + [xs(,-T),[x,]]) 



dTD{T){[x,xs{-r)l}. 



(4-6) 



Thus, we expect {A)t represents the quantum counter-part of the phase volume 
contraction. 

In the steady state limit, we should take the limit of t ^ oo as 

;>oo 

(p)ss = lim {p)t = / dttIs[pcc^epH{t)]. (4-7) 
The contribution of the first term on the right hand side of (|3-15p is 

(P)ss = — r dt r dX{ps{-ihX)pH{t))o<„ (4-8) 
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which is the conventional Green-Kubo formula. The contribution of the second term 
on the right hand side of (j3-15p is 

poo p{5 

mi = / dt dX{{xs{-ihX),ps{-ihX)}pH{t)),^, (4-9) 

Jo Jo 

where we have used (p(t))eq = 0. The contribution of the third term on the right 
hand side of (|3-15|) is 

^)SS=^ P dt P dT{b^{T) d\{ps{-ih\)xs{-f)pH{t))e^ 

-D+{t) / dX{ps{-ihX)xs{-T)pH{t))eci}. (4-10) 
Jo 

The contribution of the fourth term on the right hand side of (|3-15p is 

(P)SS = ^J y [D+{T){[x,Xs{-T)]pH{t))ec^ " D-{t) {[x , X s {-f)]pH {t)U) 

(4-11) 

Thus, we obtain the generalized Green formula 

4 

{Phs = (4-12) 

i=l 

It is obvious that three terms ^i=2(p)ss represents the nonlinear correction to the 
conventional Green-Kubo formula, in which the effect of the external force appears 

(2) (3) 

through the time evolution of puit)- We also note that (p)ss and (p)gg essentially 
come from non-dissipative parts, though pnit) should involve dissipative effects. So 
far, there is no approximation once we start from the basic equations presented in 
section 2. 

It should be noted that the current pnit) in Eqs. (j4-8p - ()4-lip can be replaced 
by Apnit) = pnit) — Ph{oo) if ph{oo) is finite because of Eq. (j3T6p . 

§5. Simple example 

In the previous section, we have present formal representations of generalized 
Green-Kubo formulae, but such formal expressions might be insufficient to demon- 
strate its relevancy. In this section, we demonstrate what the result is in the case 
of V{x) = 0. It should be noted that the model is exactly solvable for the harmonic 
potential but such a case there is no steady current of the particle because of the 
trap of the particle in the potential. 

If we assume V{x) = 0, xs{—t), xs{—ifiX) and ps{—ifiX) are respectively written 

as 

p 

xs{-t)=x-\ T, (5-1) 

m 

p 

xsi—ihX) = X — ihX — , (5'2) 
m 

Ps{-ihX) =p. (5-3) 
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From Eq. (j5-3p we immediately obtain 



d\ps{-ihX) = Pp. (5-4) 







With the aid pf Eqs.dMD and (iBaTj) . thus, we obtain 



(P)Si = ^ dt dX{p%,e-^^' = ^, (5-5) 
JO Jo ^7 



where we have used (p^)eq = mkT. This is the result from the conventional Green- 
Kubo formula. We also need to stress that this result is identical to the exact solution 
of pnit oo) in Eq. ()B-lip without any statistical average. Thus, we expect that 
the contributions from Eqs. (|4-9p - ()4-lip are zero in this simple example. Indeed, it 
is quite easy to prove the above statement. 

Let us evaluate the contribution of Eq. ()4-9p . From (j5-2p and (j5-3p , we obtain 



ihX 

{xs(—ifiX),ps(—ifiX)} = xp + px — 2 p^. (5-6) 

m 

Therefore, we immediately obtain 

{{xs{-ihX),ps{-ihX)}pH{t))cq = 0, (5-7) 

where we have used {{xs{—T),ps{—ihX)})eq = and function containing odd powers 
of X or p becomes traceless. Here we note pnit) is given by Eq. ()B-lip . 

Let us evaluate the integral terms on the right hand side of (|2-16p . By using 
(|5-ip we can write 

T 

[X,{xs(-T),peq}] = [x, {x, Peq}] H [x,{p,Pcq}]- (5-8) 

m 

Therefore, we directly obtain 

tr5{[x, {xs{-T),Peq}]pH{t)} = 0, (5-9) 

where we have used tr5'{[x, {xsi—T),Peq}] = and Eq. (|B-lip . 
Similarly, from 

T 

[X, [x(-r),/)cq]] = [X, [x,Peq\] H [x, [p, Pcq\] (5-10) 

m 

we immediately obtain 

tr5{[x, [xs{-T),peci]]pH{t)} = 0, (5-11) 

where we have used tr^^x, [xsi—T), Pcq\] = and Eq. (|B-lip . 

Therefore, the generalized Green-Kubo formula for the motion without potential 
is reduced to the result obtained by conventional Green-Kubo formula as 

(P)ss = (5-12) 
27 

To know the nonlinear contributions in Eqs. (|4-9p - (j4-lip explicitly, we need to in- 
troduce nonlinear effects of potential or the interaction between particles. 
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§6. Discussion and Conclusion 



We have obtained the generahzed Green-Kubo formula. The final expression 
should be nearly equal to (j4-8p - (j4-lip . We also verify the validity of our formulation 
in the simplest case for a free Brownian particle under the external force F^^, where 
the result is obtained from the conventional Green-Kubo formula. 

There are couple of unsolved questions to be answered, (i) The generalized 
Green-Kubo formula is not directly related to that we have obtained for classical 
systems^ where is sum of P{Hs)eq — {^)eq! where (yl)eq is the classical phase 
volume contraction, (ii) SaitcP^ found that {Hs) (more precisely, quantum version 
of Rayleigh's dissipation function) is directly related to quantum version of Harada- 
Sasa relation.'^ However, the connection between my formulation and quantum 
Harada-Sasa relation is not clear, (iii) We believe that it is straightforward to derive 
the integral fluctuation theorem in this context without using time-reversed path. 
This is the next task, (iv) Closely related to the above, how to understand generalized 
Onsager-Casmir relation in this context 'i® (v) How to apply this formulation to the 
case of microscopic time irreversible quantum systems'i'13 (vi) We have analyzed 
a case of steady external force, but we should extend the formulation for the case 
of time-dependent external field. This is indeed necessary to discuss fluctuation- 
dissipation relation and its violation, (vii) We should analyze the case of nonlinear 
potentials to clarify the correction of conventional Green-Kubo formula. In this case, 
we cannot obtain the exact solution. 
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[A,e-^"^]=e~^"^ dXe^"^[Hs,A]e-^"^ = -ihe~^"^ dXAs{-ihX), (A-1) 



where A is any observable. One can generalize this identity to Fermionic commuta- 
tion relation or the case including double commutators as 

{A,e-^"s^ = 2e-^"sA + e-^"s / d\e^"s[Hs, A]e-^^s , (A-2) 

Jo 
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Appendix A 

Generalized Kubo's identity 



Kubo used an important identity in his paper.'^ This identity can be written as 
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[A,e-^"^Bsi-T)] = e-'^^^[A,Bs{-T)] + e'^^^ / dXe^^^[Hs,A]e-^"^ Bsi-r) 

Jo 

= e-^^^[A,Bs{-T)]-ihe-^^^ dXAs{-ihX)Bs{-T). (A-3) 

JO 

The derivation of these identities are straightfoward. First, we derive eq. ()A-2p . 
It is easy to confirm 

(A-4) 

On the other hand, it is easy to confirm the identity 

d rf^ 

■i- dXe^"^-[Hs,A]e-^"s =e^^^[Hs,A]e-^"s. (A-5) 
dp Jo 

Thus, we readily obtain eq. ()A-2j) where we have used e^^^ {A, e~^^^} 2A as 
/3 — >■ 0. The derivation of eq. ()A-3P is almost identical to the above. If we can use 

{d/dp)e(^"s[A,e-^"s Bs{-t)] = e^^^lHs, A]e-P"s Bs{-t) ande^"s[A,e-^"s Bs{-t)] 
[A^Bsi—r)] as /? ^ 0, we readily obtain eq. ()A-3p . 

Appendix B 

Explicit calculation ofpnit) 



It is possible to obtain the exact solution of pnit) if the potential is harmonic 
(??). The formal solution pnit) is written as 



^ t 



PH{t) = e'^'p = Y,-,my, (B-i) 



nl 

n=0 



where iC is given by (j4-4p . Thus, the most important process to obtain pnit) is to 
obtain iCp which consists of the five terms as iC = Ylf=i The first of iCp is 



The second term of iCp is 
The third term is given by 



iCiP=^[Hs,p]=0. (B-2) 



iC2P=^[Hc,p] = -2ICx. (B-3) 



iCsp = -'-^[x,p]=Fe.. (B-4) 
The fourth term is the most complicated, which is given by 

iCiP = dTD{T){xs{-T), [x,p]} = -^J^ dTD{T)xs{-T) 

= 1 r dTD{T)-— r dTD{T)T 

= -2jp + 2JCx, (B-5) 
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where we have used (15-ip for the second equality, and 



/ cItD^t) = 2h hm / dcu dTj{uj)e 
Jo ^-"^Jo Jo 



= 2h\\m I duj^^^^^ = 2hK. (B-6) 

We should emphasize that the contribution of counter Hamilitonian ()B-3P is cancelled 
from the contribution of Eq. (|B-6p . The five contribution to iCp is given by 

1 

iC-^P ^ dTDi{T)[xs{-T), [x,p]] = 0. (B-7) 

Thus, from Eqs. (|B-2p - ()B-7p we obtain 

iCp = -2jp + Fe., (B-8) 

or equivalently 

iCAp = -2jAp, (B-9) 

where 

Ap = p-^. (B-10) 

Thus, we obtain (i£)"Z\p = (-27)"Z\p. 
This leads to 

PHit) = t'^my = If + - If ) (i^-ii) 
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